Percolation in Networks with Voids and Bottlenecks 
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A general method is proposed for predicting the asymptotic percolation threshold of networks 
with bottlenecks, in the limit that the sub-net mesh size goes to zero. The validity of this method 
is tested for bond percolation on filled checkerboard and "stack-of-triangle" lattices. Thresholds for 
the checkerboard lattices of different mesh sizes are estimated using the gradient percolation method, 
while for the triangular system they are found exactly using the triangle-triangle transformation. 
The values of the thresholds approach the asymptotic values of 0.64222 and 0.53993 respectively 
as the mesh is made finer, consistent with a direct determination based upon the predicted critical 
corner-connection probability. 

PACS numbers: 64.60.ah, 64.60.De, 05.50.+q 



I. INTRODUCTION 

Percolation concerns the formation of long-range con- 
nectivity in random systems [jj. It has a wide range 
of application in problems in physics and engineering, 
including such topics as conductivity and magnetism in 
random systems, fluid flow in porous media 0], epidemics 
and clusters in complex networks Q, analysis of water 
structure and gelation in polymer systems @. To 
study this phenomenon, one typically models the net- 
work by a regular lattice made random by independently 
making sites or bonds occupied with probability p. At 
a critical threshold p c , for a given lattice and percola- 
tion type (site, bond), percolation takes place. Find- 
ing that threshold exactly or numerically to high pre- 
cision is essential to studying the percolation problem 
on a particular lattice, and has been the subject of nu- 
merous works over the years (recent works include Refs. 

jl| l^i^'Q^ IS IS IS HE UE IS IS El HE HH, 

In this paper we investigate the percolation character- 
istics of networks with bottlenecks. That is, we consider 
models in which we increase the number of internal bonds 
within a sub-net while keeping the number of contact 
points between sub-nets constant. We want to find how 
p c depends upon the mesh size in the sub-nets and in 
particular how it behaves as the mesh size goes to zero. 
Studying such systems should give insight on the behav- 
ior of real systems with bottlenecks, like traffic networks, 
electric power transmission networks, and ecological sys- 
tems. It is also interesting from a theoretical point of 
view because it interrelates the percolation characteris- 
tics of the sub-net and the entire network. 

An interesting class of such systems includes lattices 
with an ordered series of vacated areas within them. Ex- 
amples include the filled checkerboard lattices (Fig. [1]) 
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and the "stack-of-triangles" (Fig. [2]) . The latter can be 
built by partitioning the triangular lattice into triangu- 
lar blocks of dimension L, and alternately vacating those 
blocks. These internal blocks of length L correspond 
to the sub-nets, which contact other sub-nets through 
the three contact points at their corners. The checker- 
board lattice is the square-lattice analog of the stack-of- 
triangles lattice, where sub-nets are L x L square lattices 
which contact the other sub-nets via four contact points. 
Note, for the stack-of-triangles sub-nets, we also use the 
L x L designation, here to indicate L bonds on the base 
and the sides. 

The problem of finding the bond percolation thresh- 
old can be solved exactly for the stack-of-triangles lat- 
tice because it fits into a class of self-dual arrangements 
of triangles, and the triangle-triangle transformation (a 
generalization of the star-triangle transformation) can be 
used to write down equations for its percolation thresh- 
old [13, Hl| . This approach leads to an algebraic equa- 
tion which can be solved using numerical root-finding 
methods. However due to lack of self-duality in the filled 
checkerboard lattices, no exact solution can be obtained 
for their thresholds. 

It is of interest and of practical importance to inves- 
tigate the limiting behavior of systems with sub-nets of 
an infinite number of bonds, i.e., systems where the size 
of sub-nets is orders of magnitude larger than the size of 
a single bond in the system, or equivalently, where the 
mesh size of the lattice compared to the sub-net size be- 
comes small. Due to reduced connectivity, these systems 
will percolate at a higher occupation probability than a 
similar regular lattice. The limiting percolation thresh- 
old for infinite sub-nets is counter-intuitively non-unity, 
and is argued to be governed by the connectedness of con- 
tact points to the infinite percolating clusters within sub- 
nets. This argument leads to a simple criterion linking 
the threshold to the probability that the corners connect 
to the giant cluster in the center of the sub-net. 

In this work, the limiting threshold value is computed 
for bond percolation on the stack-of-triangles and filled 
checkerboard lattices using this new criterion. Percola- 
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FIG. 1: Checkerboard lattices with sub-nets of finite sizes. 
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FIG. 2: Stack-of-triangles lattice with sub-nets of finite sizes. 



tion thresholds are also found for a series of lattices of 
finite sub-net sizes. For the stack-of-triangles lattices, 
most percolation thresholds are evaluated analytically 
using the triangle-triangle transformation method, while 
for filled checkerboard lattices, the gradient percolation 
method [2!| is used. The limiting values of 0.53993 and 
0.64222 are found for percolation thresholds of stack-of- 
triangles and checkerboard lattices respectively, which 
are both in good agreement with the values extrapolated 
for the corresponding lattices of finite sub-net sizes. 

We note that there are some similarities between this 
work and studies done on the fractal Sicrpihski gaskets 
(triangular) [3(| and carpets (square), but in the case 
of the Sierpihski models, the sub- nets arc repeated in 
a hierarchical fashion while here they are not. For the 
Sierpihski gasket, which is effectively all corners, the per- 
colation threshold is known to be 1 (3lj . For Sierpihski 
gaskets of a finite number of generations, the formulae for 
the corner connectivities can be found exactly through 
recursion [32j, while here they cannot. Recently another 
hierarchical model with bottlenecks, the so-called Apol- 
lonian networks, which are related to duals of Sierpinski 
networks, has also been introduced [331 ]. In this model, 
the percolation threshold goes to zero as the system size 
goes to infinity. 



II. THEORY 

Let p be the probability that a bond in the system is 
occupied. Consider a network with sub-nets of infinitely 
fine mesh, each individually percolating (in the sense of 
forming "infinite" clusters but not necessarily connecting 
the corners) at p c , Sl and denote the overall bond perco- 
lation threshold of the entire network to be p c . n . It is 
obvious that p c , s < p c ,n, due to reduced connectivity in 
the entire network compared to connectivity in individ- 
ual sub- nets. For p c s < p < p CjTl , an infinite cluster will 
form within each sub-net with probability 1. However, 
the entire network will not percolate, because a suffi- 
cient number of connections has not yet been established 
between the contact points at the corners and central 
infinite clusters. 

Now we construct an auxiliary lattice by connecting 
the contact points to the center of each subnet, which 
represents the central infinite cluster contracted into a 
single site. The occupation probability of a bond on this 
auxiliary lattice is the probability that the contact point 
is connected to the central infinite cluster of the sub-net. 
Percolation of this auxiliary lattice is equivalent to the 
percolation of the entire network. That is, if this auxil- 
iary lattice percolates at a threshold p c . a , the percolation 
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FIG. 3: (Color online.) Stack-of-triangles lattice and its aux- 
iliary lattice. The filled blue (dark) triangles represent the 
sub-net, and the yellow honeycomb lattice represents the ef- 
fective auxiliary lattice. 



FIG. 4: (Color online.) Auxiliary lattice of the checkerboard 
lattice. The blue (dark) colored areas represent the subnets, 
and the double-bond square lattice (diagonals) represents the 
auxiliary lattice. 



threshold of the entire network will be determined by: 



oo, corner 



(Pc 



(1) 



where P oc _ corner (p) gives the probability that the corner 
of the sub-net is connected to the central infinite cluster 
given that the single occupation probability is p. 

In general no analytical expression exists for 
-Poo,comer(p), even for simple lattices such as the trian- 



gular and square lattices, and P r 
uated by simulation. 



oo, corner 



(p) must be eval- 



A. Stack-of- Triangles Lattice 



B. Filled Checkerboard Lattice 

Unlike the stack-of-triangles lattice, there is no exact 
solution for percolation threshold of the checkerboard 
lattice for finite sub-nets because no duality argument 
can be made for such lattices. However once again an 
auxiliary lattice approach can be used to find a crite- 
rion for the asymptotic value of percolation threshold. 
Fig. [4] depicts the corresponding auxiliary lattice for a 
checkerboard lattice, which is simply the square lattice 
with double bonds in series. This lattice percolates at 
p CiCl = 1/V5« 0.707107. Thus for the infinite sub-net 
p Cj „ will be determined by: 



Fig. [3] shows a limiting stack-of-triangles lattice 
where each shaded triangle represents a sub-net of 
infinitely many bonds. The contact points are the 
corners of the triangular sub-nets. As shown in 
Fig. [3l the auxiliary lattice of the stack-of-triangles 
lattice is the honeycomb lattice, which percolates at 



p Cia = 1 - 2sin(7r/18) « 0.652704 [3J. Thus the asymp- 
totic percolation threshold p c 
will be determined by: 



of the stack-of-triangles 



oo, corner (Pc,n) a/2 



(4) 



It is interesting to note that there exists another regular 
lattice — the "martini" lattice — for which the bond 
threshold is also exactly 1/V2 @. However, that lattice 
does not appear to relate to a network construction as 
the double-square lattice does. 



P, 



oo, corner 



(p c ,n) = 1 - 2 sin 



(2) 



Because the stack-of-triangles lattice is made up of tri- 
angular cells in a self-dual arrangement, its percolation 
threshold can be found exactly using the triangle-triangle 
transformation [27], [28[. Denoting the corners of a sin- 
gle triangular sub-net with A, B and C, the percolation 
threshold of the entire lattice is determined by the solu- 
tion of the following equation: 



P(ABC) = P(ABC) 



(3) 



where P(ABC) is the probability that A, B and C are 
all connected, and P(ABC) is the probability that none 
of them are connected. Eq. (J3j) gives rise to an algebraic 
equation which can be solved for the exact percolation 
threshold of the lattices of different sub- net sizes. 



III. METHODS 



A. Percolation Threshold of Systems of finite-sized 
sub-nets 



For the checkerboard lattice, we estimate the bond 
percolation thresholds using the gradient percolation 
method [29j |. In this method, a gradient of occupation 
probability is applied to the lattice, such that bonds are 
occupied according to the local probability determined 
by this gradient. A self-avoiding hull-generating walk is 
then made on the lattice according to the rule that an 
occupied bond will reflect the walk while a vacant bond 
will be traversed by the walk. For a finite gradient, this 
walk can be continued infinitely by replicating the origi- 
nal lattice in the direction perpendicular to the gradient 
using periodic boundary conditions. Such a walk will 



map out the boundary between the percolating and non- 
percolating regions, and the average value of occupation 
probability during the walk will be a measure of the per- 
colation threshold. Because all bonds are occupied or 
vacated independent of each other, this average proba- 
bility can be estimated as (35l |: 



It is particularly straightforward to implement this algo- 
rithm to bond percolation on a square lattice, and the 
checkerboard lattice can be simulated by making some 
of the square-lattice bonds permanently vacant. Walks 
are carried out in a horizontal-vertical direction and the 
original lattice is rotated 45°. 

We applied this approach to checkerboard lattices of 
different block sizes. Fig. \5\ and Fig. [H] show the corre- 
sponding setups for lattices with 2x2 and 4x4 vacancies, 
where the lattice bonds are represented as dashed diag- 
onal lines and solid horizontal and vertical lines show 
where the walk goes. Circles indicate the centers of per- 
manently vacant bonds. It should be emphasized that 
permanently vacated bonds are not counted in Eq. ([5]) 
even if they are visited by the walk. 

The percolation threshold of stack-of-triangles lattices 
of finite sub- net size were calculated using Eq. ([3]). If the 
occupation probability is p and q = 1 —p, one can express 
P(ABC) and P(ABC) as: 

3n(n+l)/2 

P(ABC) = <HMpY™ (n+1)/2 ~ l (6) 

3n(ra+l)/2 

P(ABC) = ^(n>jy? 3 ' l(,!+1)/2 " i (7) 

i=0 

where n denotes the number of bonds per side of the sub- 
net, 4>{n, i) denotes the number of configurations of an nx 
n triangular block with precisely i occupied bonds where 
the A, B and C are connected to each other and 
denotes the number of configurations where none of these 
points are connected. There appears to be no closed- 
form combinatorial expression for i) and ip(n, i), and 
we determined them by exhaustive search of all possible 
configurations. 

B. Estimation Of Poo, corner 

As mentioned in Section [H] the asymptotic value 
of percolation threshold p c ^ n can be calculated using 
Eq. ([4]). However there is no analytical expression for 
-Poo,comer(.p) 5 hence it must be characterized by simula- 
tion. In order to do that, the size distribution of clusters 
connected to the corner must be found for different val- 
ues of p > Pes- Cluster sizes are defined in terms of the 
number of sites in the cluster. 

In order to isolate the cluster connected to the cor- 
ner, a first-in- first-out Leath or growth algorithm is used 
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FIG. 5: Representation of checkerboard lattices for simula- 
tion with gradient method. The original bond lattice is rep- 
resented by dashed diagonal lines, while lattice on which the 
walk goes is vertical and horizontal. Open circles mark bonds 
that are permanently vacant. 
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FIG. 6: Checkerboard lattice with 4x4 vacancies, with de- 
scription the same as in Fig. [5] 



starting from the corner. In FIFO algorithm, the neigh- 
bors of every unvisited site are investigated before going 
to neighbors of the neighbors, so that clusters grow in 
a circular front. Compared to last-in-first-out algorithm 
used in recursive programming, this algorithm performs 
better for p>p c , s because it explores the space in a more 
compact way. 

At each run, the size of the cluster connected to the cor- 
ner is evaluated using the FIFO growth algorithm. In or- 
der to get better statistics, clusters with sizes between 2 1 
and 2 l+1 — 1 are counted to be in the i— th bin. Because 
simulations are always run on a finite system, there is an 
ambiguity on how to define the infinite cluster. However, 
when p>p c .s, the infinite cluster occupies almost the en- 
tire lattice, and the finite-size clusters are quite small on 
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TABLE I: Percolation threshold for checkerboard lattices of 
different sub-net sizes: "Exact result, ''from gradient percola- 
tion simulations, c from corner simulations using Eq. Q. 



B. 



Percolation Threshold of The Stack-of-triangles 
Lattice 



Sub- net size 


Estimated p c ,„ 


lxl 


0.5° 


2x2 


0.596303 ± 0.000001 b 


4x4 


0.633685 ± O.OOOOOQ 6 


8x8 


0.642318 ± 0.000005 6 


16 x 16 


0.64237 ± 0.00001'' 


32 x 32 


0.64219 ± 0.00002 6 


oo 


0.642216 ± 0.00001 c 



average. This effect becomes more and more profound 
as p increases, and the expected number of clusters in a 
specific bin becomes smaller and smaller. Consequently, 
larger bins will effectively contain no clusters, except the 
bin corresponding to cluster sizes comparable to the size 
of the entire system. Thus there is no need to set a cutoff 
value for defining an infinite cluster. Fig. [7] depicts the 
size distribution of clusters connected to the corner ob- 
tained for 1024 x 1024 triangular lattice at (a): p — 0.40 
and (b): p = 0.55 after 10 4 independent runs. As it is ob- 
served, there is a clear gap between bins corresponding to 
small clusters and the bin corresponding to the spanning 
infinite cluster even for small values of p, which clearly 
demonstrates that the largest nonempty bin corresponds 
to infinite percolating clusters connected to the corner. 
The fraction of such clusters connected to the corner is 
an estimate of Poo, comer (p)- 

In the simulations, we used the four-offset shift-register 
random-number generator R(471, 1586, 6988, 9689) de- 
scribed in Ref. [361 ] . 



IV. RESULTS AND DISCUSSION 

A. Gradient Percolation Data 

The gradient percolation method was used to estimate 
the bond percolation threshold of checkerboard lattices of 
five different sub-net sizes, i.e., 2 x 2, 4 x 4, 8 x 8, 16 x 16 
and 32 x 32. For each lattice, six values of the gradi- 
ent were used, and simulations were run for 10 10 to 10 12 
steps for each gradient value in order to assure that the 
estimated percolation thresholds are accurate to at least 
five significant digits. The gradient was applied at an 



angle of 45° relative to the original lattice. Figures 8(a) 
8(c) depict typical simulation results. Measured perco- 
lation thresholds for finite gradients were extrapolated 
to estimate the percolation threshold as L — > oo. Our 
simulations show that p c fits fairly linearly when plotted 
against 1/L. Table U gives these estimated percolation 
thresholds. 



As mentioned in Sections III Al and IIII A[ the percola- 
tion threshold of stack-of-triangles lattice can be deter- 
mined by Eq. ([3]) . Table [TT| summarizes the correspond- 
ing polynomial expressions and their relevant roots for 
lattices having 1, 2, 3 and 4 triangles per edge. These 
polynomials give the </>(«., i) and ip(n, i) in Eqs. ^ and 
© for n = 1,2,3, 4 and i = 0,1,.. .,3n(n + l)/2. We 
showpo = P(ABC), pi = P(ABC) (the probability that 
a given pair of vertices are connected together and not 
connected to the third vertex), andp3 = P(ABC). These 
quantities satisfy p + 3p 2 +Pz = 1. Then we use Eq. © 
to solve for p Cj „ numerically. 

We also show in Table UT1 the values of po, Pi and P3 
evaluated at the p c , n - Interestingly, as n increases, po 
at first increases somewhat but then tends back to its 
original value at n — 1, reflecting the fact that the con- 
nectivity of the infinitely fine mesh triangle is identical 
to that of the critical honeycomb lattice, which is iden- 
tical to the connectivity of the simple triangular lattice 
according to the usual star-triangle arguments. 

It is not possible to perform this exact enumeration for 
larger sub-nets, so we used gradient percolation method 
to evaluate p c for 5x5. (To create the triangular bond 
system on a square bond lattice, alternating horizon- 
tal bonds are made permanently occupied.) The final 
threshold results are summarized in Table [TTTT 



C. Estimation Of Poo,corncr(p) 

1. Square Lattice 

The cluster growth algorithm was used to estimate 
-Poo, comer (p) for different values of p. Simulations were 
run on a 2048 x 2048 square lattice. For each value 
of p > 1/2, 10 5 independent runs were performed and 
Poo, comer was estimated by considering the fraction of 
clusters falling into the largest nonempty bin as described 
in Section UlIBI Fig. [5] demonstrates the resulting curve 
for the square lattice. In order to solve Eq. (01, a cubic 
spline with natural boundary conditions was used for in- 
terpolation, and an initial estimate of p CjTl was obtained 
to be 0.6432. The standard deviation of Poo, comer (p) 



scales as 0(1/ y/N) where N is the number of indepen- 
dent simulation used for its estimation, so that N = 10 5 
will give us an accuracy in Poo,corner(p) of about two sig- 
nificant figures. 

In order to increase the accuracy in our estimate, 
further simulations were performed at the vicinity of 
p = 0.6432 for N = 10 10 trials with a lower cut-off size, 
and p Cj „ was found to be 0.642216±0.00001. This number 
is in good agreement with percolation thresholds given in 
Table HI Note that p c n for the 16 x 16 sub-net checker- 
board lattice actually overshoots the value 0.642216 for 
the infinite sub-net and then drops to the final value. 
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FIG. 7: Cluster size distribution for clusters originating from the corner of triangular lattice at (a): p = 0.40 and (b): p = 0.55, 
after 10 4 independent runs. 




0.643 

■a 



o 0.6425 




FIG. 8: Gradient percolation data for checkerboard lattices with: (a) 2x2, (b) 4x4, and (c) 8 x 8 sub-net sizes. Squares 
correspond to the bond percolation thresholds estimated for different values of probability gradients; error bars are smaller 
than the symbol size. 



This non-monotonic behavior is surprising at first and ters !l] n s ~ aexp(— bs 1 / 2 ) with In a = —7.0429 and 
presumably is due to some interplay between the vari- b = —0.8177. 
ous corner connection probabilities that occurs for finite 
system. 

At the threshold p Cj „ = 0.642216, we found that the 
number of corner clusters containing s sites for large 
s behaves in the expected way for supercritical clus- 
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TABLE II: Exact enumeration polynomials and p c — p Cj „ for the stack-of-triangles lattices 



Sub-net 


lxl (simple triangular lattice) 2 x 2 (3 up triangles or 9 bonds per sub-net) 


p 3 = P(ABC) 
P2=P (ABC) 
po = P(ABC) 

Pc 

Po(Pc) =P3(Pc) 
P2(Pc) 


p 3 + 3p 2 g 9p 4 g b + 57p 5 g 4 + 63pV + 33p Y g 2 + 9p 8 g + p H 
pq 2 p 2 q 7 + 10p 3 g 6 + 32pV + 22p 5 g 4 + 7p 6 g 3 + pV 
(p + g) 3 - p 3 - 3p 2 (p + g) 9 -pa - 3p 2 
0.34729635533 0.47162878827 
0.27806614328 0.28488908000 
0.14795590448 0.14340728000 


Sub-net 


3 x 3 (6 "up" triangles or 18 bonds per sub-net) 


p 3 = P(ABC) 
p 2 = P (ABC) 
po = P(ABC) 

Pc 

Po(Pc) = P3(pc) 
P2(Pc) 


29p B g i2 +468p'g ii + 3015p 8 g iu + 9648p u g a + 16119p lu g 8 + 17076p il g' + 12638p i: V 

+6810p 13 g 5 + 2694p 14 g 4 + 768p 15 g 3 + 150p 16 g 2 + 18p 17 g + p 18 

p 3 g 15 + 21p 4 g 14 + 202p 5 g 13 + 1125p 6 g 12 + 3840p 7 g 11 + 7956p 8 g 10 + 9697p 9 g 9 

+7821p 10 g 8 + 4484p n g 7 + 1879p 12 g 6 + 572p 13 g 5 + 121p 14 g 4 + 16p 15 g 3 +p 16 g 2 

(p + g) 18 -p3-3p 2 

0.50907779266 

0.28322276251 

0.14451815833 


Sub-net 


4x4 (10 "up" triangles or 30 bonds per sub-net) 


p 3 = P(ABC) 
P2 = P (ABC) 
po = P(ABC) 

Pc 

Po(Pc) = P'i{Pc) 
P2(Pc) 


99p 8 g 22 + 2900p 9 g 21 + 38535p iu g 2U + 305436p ii g 19 + 1598501p i2 g 18 

+5790150p 13 g 17 + 14901222p 14 g 16 + 27985060p 15 g 15 + 39969432p 16 g 14 

+45060150p 17 g 13 + 41218818p 18 g 12 + 31162896p 19 g 11 + 19685874p 20 g 10 

+10440740p 21 g 9 + 4647369p 22 g 8 + 1727208p 23 g 7 + 530552p 24 g 6 + 132528p 25 g 5 

+26265p 26 g 4 + 3976p 27 g 3 + 432p 28 g 2 + 30p 29 g +p 30 

p 4 g 26 + 36p 5 g 25 + 613p 6 g 24 + 65 33p 7 g 23 + 486 43p 8 g 22 + 26 7261p 9 g 21 

+1114020p 10 g 20 + 3563 824p 11 g 19 + 8766414p 12 g 18 + 16564475p 13 g 17 

+24187447p 14 g 16 + 27879685p 15 g 15 + 25987202p 16 g 14 + 19980934p 17 g 13 

+12843832p 18 g 12 +6950714p 19 g n + 3170022p 20 g 10 + 1212944p 21 g 9 

+385509p 22 g 8 + 100140p 23 g 7 + 20744p 24 g 6 + 3300p 25 g 5 + 379p 26 g 4 + 28p 27 g 3 + p 28 g 2 

(p + q) 30 -P3 -3p 2 

0.52436482243 
0.28153957013 
0.14564028658 



TABLE III: Percolation threshold for stack-of-triangles lat- 
tices of different sub-net sizes: "From Eq. ((3} using exact 
expressions for po and P3 from Table ITU "from gradient simu- 
lation, c Corner simulation using Eq. ([2]). 



Sub-net size 


Estimated p c 


lxl 


0.347296355° 


2x2 


0.471628788° 


3x3 


0.509077793° 


4x4 


0.524364822° 


5x5 


0.5315976 ± 0.000001 6 


oo 


0.53993 ± 0.00001 c 




0.6 0.7 0.8 0.9 1 

Occupation Probability 



FIG. 9: P c 



oc, corner 



(p) for the square lattice. 



2. Triangular lattice 



The cluster growth algorithm was applied to find the 
size distribution of clusters connected to the corner of 
a 1024 x 1024 triangular lattice. For each value of p, 
10 4 independent runs were performed and Poo, comer (p) 
was evaluated. Fig. [10] depicts the results. The root 
of Eq. ([J) was determined by cubic spline to be around 



0.539. Further simulations were performed around this 
value with TV = 10 10 runs for each p, yielding p Cj „ = 
0.539933 ± 0.00001. This value is also in good agreement 
with values given in Table [ill and shows fast convergence 
as sub-net size increases. 
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o.2- y 



0.4 0.5 0.6 0.7 0.8 0.9 1 

Occupation Probability 



FIG. 10: Poo, comer (p) for the triangular lattice. 




FIG. 11: A particular realization of the 8x8 checkerboard 
lattice at its critical point p c ,„ = 0.642318 (as given in Table 
HJ, showing in black the sites wetted by the largest cluster. 
The total lattice is 128 x 128 sites, with periodic b.c. 



V. DISCUSSION 

We have shown that the percolation threshold of 
checkerboard and stack-of-triangle systems approach val- 
ues less than 1 as the mesh spacing in the sub-nets goes 
to zero. In that limit, the threshold can be found by 
finding the value of p such that the probability a corner 
vertex is connected to the infinite cluster Poo, comer equals 
1/V2 and 1 — 2sin(7r/18), respectively, based upon the 
equivalence with the double-bond square and bond hon- 
eycomb lattices. The main results of our analysis and 
simulations are summarized in Tables [J and IIIII 

For the case of the checkerboard, we notice a rather 
interesting and unexpected situation in which the thresh- 
old p c ,n slightly overshoots the infinite-sub-net value and 
then decreases as the mesh size increases. The thresh- 
old here is governed by a complicated interplay of con- 
nection probabilities for each square, and evidently for 
intermediate sized systems it is somewhat harder to con- 
nect the corners than for larger ones, and this leads to 




FIG. 12: Average density of the infinite (largest) cluster in 
a 64 x 64 square, at the checkerboard critical point p c ,„ = 
0.642216. At the corners the probability drops to l/y/2 = 
0.707 according to Eq. Q. 



a larger threshold. In the case of the triangular lattice, 
where there are fewer connection configurations between 
the three vertices of one triangle (namely, just po, pi and 
P3), the value of p c _ n appears to grow monotonically. 

To illustrate the general behavior of the systems, we 
show a typical critical cluster for the 8x8 checkerboard 
system in Fig. QT] It can be seen that the checkerboard 
squares the cluster touches are mostly filled, since the 
threshold p c>n — 0.642318 is so much larger than the 
square lattice's threshold p CtS = 0.5. 

In Fig. [T2] we show the average density of "infinite" 
(large) clusters in a single 64 x 64 square at the checker- 
board criticality of p c , n = 0.642216, in which case the 
density drops to l/v2 at the corners. In Fig.[l3]we show 
the corresponding densities conditional on the require- 
ment that the cluster simultaneously touches all four cor- 
ners, so that the density now goes to 1 at the corners and 
drops to a somewhat lower value in the center because 
not every site in the system belongs to the spanning clus- 
ter. Similar plots can be made of clusters touching 1, 2, 
or 3 corners. At the sub-net critical point p c .s, the first 
two cases can be solved exactly and satisfy a factoriza- 
tion condition [37], [38| , but this result does not apply at 
the higher p c , n . 

The ideas discussed in this paper apply to any sys- 
tem with regular bottlenecks. Another example is the 
kagome lattice with the triangles filled with a finer-mesh 
triangular lattice; this system is studied in Ref. [39| . 
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FIG. 13: Average density of clusters that simultaneously 
touch all four corners at the checkerboard criticality p c , n = 
0.642216 for a single 64 x 64 square. 
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